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Abstract 

At proposed high energy linear e + e~ colliders a large number of W bosons and top 
quarks will be produced. We evaluate the potential precision to which the decay branching 
ratios into the various quark species can be measured, implying also the determination 
of the respective CKM matrix elements. Crucial is the identification of the individual 
quark flavours, which can be achieved independent of QCD models. For transitions 
involving up quarks the accuracy is of the same order of magnitude as has been reached 
in hadron decays. We estimate that for charm transitions a precision can be reached 
that is superior to current and projected traditional kinds of measurements. The t — > b 
determination will be significantly improved, and for the first time a direct measurement 
of the t — > s transition can be made. In all cases such a determination is complementary 
to the traditional way of extracting the CKM matrix elements. 



1 Introduction 



There are nineteen fundamental parameters which cannot be derived from first principles 
in the Standard Model [I],|2| with three fermion generations and massless neutrinos. Three 
of these are the mixing angles of the Cabibbo-Kobayashi-Maskawa (CKM) matrix || 
which express the probability for a transition between different quark species through 
charged weak currents. As yet it is unclear how the number and existence of generations 
can be derived from first principles of an underlying theory. Thus, no prescription exists 
on how to calculate the values of the CKM matrix elements. Still, it is important to 
measure their values to the best precision possible. Up to now the CKM matrix elements 
have been derived from hadron decays using QCD symmetries to extend the theoretical 
analysis into a low energy scale. The precision of such measurements is therefore often 
limited by uncertainties reflecting theoretical models and the assumptions invoked. In 
particular, the transitions involving top quarks can only be indirectly determined, using 
either B^, B s mixing or b — > S7 decays. 

The direct observation of the decay products of W bosons and top quarks offers a 
complementary way to determine the CKM matrix elements. In case of W bosons one 
may determine the decay widths of W bosons into identified pairs of quarksQ. Since the 
branching ratio of the W bosons into a specific pair of quarks (q, q') is proportional to the 
square of the CKM matrix element \V qq '\, 

B(W^qq') ex |V | 2 , 

the measurement of the decay fractions allows a direct determination of the CKM matrix 
elements. Similarly one can relate the transitions of B(t — > qW) to the CKM matrix 
elements |Vt 9 |. In practice this requires that all decay modes of the W bosons and top 
quarks can be reconstructed with a good signal to background ratio, which is the case at 
e + e~ colliders. Experimentally much more complicated is the requirement that all quark 
species can be individually identified. A first step to determine the CKM matrix elements 
in W decays has been made at LEP by deriving \V CS \ from the total hadronic W branching 
ratio j|] or more directly from the measurement of the inclusive decay fraction into charm 
quarks 0,0. In both cases the other CKM matrix elements are used as a crucial input to 
the analysis when its unitarity is assumed. Our method does not rely on these assumptions 
and provides a determination of the matrix elements that is complementary to traditional 
means. 



Quark tagging for charm and bottom is based on their well understood and 
unambiguous special properties, particularly their masses and long lifetimes. The 
identification of light quarks with a precision needed for a meaningful measurement of 
the CKM matrix elements is much more involved. Several analyses have used light 
flavour tagging methods based on model assumptions such as those used in JETSET H 
and HERWIG 0] at the price of sizeable uncertainties. A method was suggested in [10| to 



1 To simplify the text, we drop the distinction between particles and antiparticles where the meaning 
it otherwise clear. 
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determine individual up, down and strange quark tagging efficiencies from Z° data, thus 
avoiding inherently ambiguous assumptions about the process of hadronisation. Such an 
analysis was recently carried out by the OPAL collaboration using 2.8 million Z° decays 
recorded at LEP pi 



The main idea for determining the fraction of light flavours is that particles with a 
large fraction x p = 2p/E cm of the momentum, p, relative to the centre-of-mass energy, 
E cm , carry information about the primary flavour 0. At the Z° the yields of single tags, 
where just one jet is tagged by a high x p particle, and the numbers of double tags, where 
both jets are tagged, can be used to determine the tagging efficiencies rf q with minimal 
reliance on assumptions about hadronisation. Here r] q is the probability that a quark q 
leads to a particle i that has the highest momentum in an event hemisphere^]. 

The efficiencies can be applied almost directly to the decays of W bosons in order to 
determine the hadronic branching fractions 

_ B(W -> qqi) 

-"-99' ~~ TD ' 

where is the inclusive branching ratio of the W to hadrons, B(W — ► hadrons). The 
numbers of single tags and double tags from W boson decays can be expressed by the 
tagging efficiencies rj l obtained from the analysis of Z° decays and the R qq '. One obtains 
an over-constrained linear equation system which can be solved for the R qq '. With the 
current samples at LEP of some 20000 W bosons per experiment, the precision of these 
decay fractions will be fairly limited. The high luminosity at proposed linear e + e~ colliders 
both for running at the Z° and at high energies above the W pair threshold offers unique 
possibilities to pursue CKM matrix measurements with substantially higher precision. A 
similar equation system can be constructed for top decays. In this note we will outline a 
strategy and estimate the potential precision. 

The outline of the paper is as follows. In Section |2| we detail the assumptions 
on accelerator and detector performance made for this study. The formalism and 
experimental aspects of the determination of light flavour tagging efficiencies at the Z° 
peak are discussed in Section |3|. The method to determine the W boson branching ratios 
is given in Section f| and their translation into the CKM matrix elements in Section |^. 
Direct measurements using top quarks are discussed in Section ^| before concluding in 
Section [7|. 



2 Assumptions on Data and Detector Performance 

For definiteness we assume in this paper the parameters of the TESLA [IB option for a 
future e + e~ linear collider. However, our proposal is in no way restricted to this option 

2 An event hemisphere is defined by the plane perpendicular to the event thrust axis and passing 
through the origin. In this analysis, we denote hemispheres as representing quark jets, since we are 
interested in studying the evolution of primary quarks into different hadron types. 
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e c 


eb 




charm quark tag 


0.60 


0.20 


0.04 


bottom quark tag 


0.02 


0.50 


0.0008 


light quark tag (W analysis) 


0.50 


0.10 


0.99 


light quark tag (top analysis) 


n.a. 


0.01 


0.50 



Table 1: Expected efficiencies at TESLA for algorithms optimised to identify charm and 
bottom quarks. Also listed are the efficiencies assumed for light quark tagging which are 
different for the W boson and top quark analyses. 



and can be pursued at any linear collider allowing for high luminosity measurements at 
the Z° peak (GigaZ) and energies above the top pair threshold. For running at the Z° 
peak we assume an instantaneous luminosity of £zo = 7 x 10 33 sec -1 cm~ 2 [14]. For a 
nominal year of running (i.e. 100 days of full efficiency) this implies 70 fb _1 of data, or 
some factor 500 more than collected by each LEP experiment over five years. We will 
therefore assume for this study a sample of 2 x 10 9 hadronic Z° decays. At centre-of-mass 
energies of 500 GeV a luminosity of Che = 3 x 10 34 sec _1 cm -2 is expected . 



The detector capabilities which are relevant for our method to determine the CKM 
matrix elements are in particular hadron and, to a lesser extent, lepton identification. In 
addition, it is crucial to identify secondary vertices to tag charm and bottom quarks with 
high efficiency and purity. 



Bottom and charm tagging at TESLA has been discussed in detail in ||16|| . From these 
considerations we assume efficiencies and purities as given in Table |1| These identification 
potentials^ are far better than what has been achieved at LEP because of the smaller 
beam pipe and the advanced micro vertex detectors foreseen at TESLA. By selecting jets 
without a prominent secondary vertex one can also increase the purity of a light flavour 
sample. The optimisation of the purity and efficiency of such tags depends on the process. 
Different working points are used for the analyses of W boson and top quark production. 

No strong emphasis has yet been placed on hadron identification for a TESLA detector. 
However, the proposed TPC as a central detector offers the possibility to determine the 
particle species by measuring the ionisation loss. Current estimates assume a dE/dx 
resolution of 4.5%. However, improvements may be possible ||18||. The momentum 



resolution will be substantially better than at LEP, helping to reconstruct invariant 
masses of resonances. On the other hand, the acceptance for long-lived particles such 
as Kg and A is reduced due to the higher momentum of these hadrons. Their large boost 
implies that a significant fraction of these will decay outside the tracking detectors. For 
definiteness we assume the same particle identification capabilities as obtained with the 



3 The study of (l6| suggests a dependence of the tagging efficiency and purity on the centre-of-mass 
energy. However, this is due to the optimisation of the tagging algorithms at Z° energies which is then 
applied to higher energies. Significant improvements can be expected if the algorithms are specifically 
optimised at high energies p7| . For simplicity we therefore assume the efficiency and purity of the 
algorithms to be independent of energy. 
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assigned 






true 
p(P) 




A(A) 




0.790 


0.062 


0.003 


0.038 


0.005 


K± 


0.146 


0.568 


0.148 


0.017 


0.026 


P(P) 


0.040 


0.246 


0.551 


0.014 


0.081 


Kg 


0.081 


0.030 


0.007 


0.691 


0.026 


A(A) 


0.047 


0.024 


0.024 


0.128 


0.696 


efficiency 


0.487 


0.441 


0.292 


0.155 


0.135 



Table 2: Fractional compositions of the identified samples (rows) in terms of the true 
tagging particle, for x p > 0.2, taken from |19|| . The rows do not add up to unity because 
of additional contributions that are not used as tags. The bottom row gives the average 
efficiency to correctly tag a hemisphere. 



OPAL detector [^| at LEP, listed in Table H, and do not attempt to calculate the potential 
small differences in particle identification efficiencies. 



3 Tagging Efficiencies from Z° Decays 



3.1 Formalism 



In this section we briefly summarise the basic formalism to determine tagging efficiencies 
from Z° data as suggested in |H| and extrapolate a recent LEP analysis [|IT] to TESLA 
conditions. As mentioned before, charm and bottom tags can be selected with high 
efficiency and purity using secondary vertex finding. The determination of the efficiency 
for light flavours is more complicated. The light flavour tags are based on particle types 
that are easy to identify, have a significant yield and which carry information about the 
primary flavour. As discussed in []20| such types are 7r =l= , K , protons, Kg, and A. We 
consider these particles as tagging particles if they have the scaled momentum in an event 
hemisphere x p = 2p/Mzo > 0.2, where the momentum, p, is determined in the Z° rest 
frame. Lifetime information yields an excellent separation of light quark jets from those 
of bottom and charm origin. Thus, we apply the high x p tag only to those jets that have 
no apparent secondary vertex, i.e. are tagged as light flavours. We assume the efficiencies 
listed in the fourth row of Table [I]. In principle jets can fulfil in parallel the requirements 
for light, charm and bottom tagging. To avoid double counting we assume the following 
priority among tags according to the achievable purity : bottom, charm and light flavour 
tags. If a jet is tagged by more than one of those only the higher priority tag is considered. 

In determining the efficiencies ?7u ds no assumption is made about the details of the 
hadronisation process like hardness, shape of the fragmentation functions, light flavour 
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composition in the hadronisation phase or the amount of resonance production. No such 
information from QCD models like JETSET || or HERWIG is needed. The only 
assumption that is invoked is that the branching ratios of the Z° into fermion pairs are as 



predicted by the Standard Model [21]. This is consistent with the high precision results 



obtained at LEP and can, at least for some flavours, be very accurately tested at GigaZ. 



As detailed in |Tl], 21J the r] q are determined by using tags in event hemispheres of 
a hadronic Z° decay. Each event is separated into two hemispheres using the plane 
perpendicular to the thrust axis containing the interaction point. In each hemisphere a 
secondary vertex is searched for and in case of a light flavour tag the highest momentum 
particle, labelled i, subject to the requirement x p > x m i n ., some minimum value. What 
can be directly observed are the number of "single-tagged hemispheres" tagged as type 
i, labelled N i: and the number of "double-tagged events" containing a tag in both 
hemispheres, labelled Nj k , where j and k are the tagging particle types. The tagging 
probability 77* is then given by 



for a number, N q , of hemispheres which originate from a quark of type q and a number, 
N q _,i, of these with tags of type i. The event counts are related to the tagging probabilities 
by: 



Ni 



2 £ rf q R q 



iV z h ad . , 

n q=d,u,s,c,b 

^d ^jpj = (2-<y Pz o ^ k q R q , (2) 

Z° g=d,u,s,c,b 

where Sjk = 1 if j = k and zero otherwise and iV z 1 o d is the number of hadronic Z° decays. 
Note that the rf q include possible distortions due to detector effects. The parameter p z o 
will not be unity if there are correlations between the tagging probabilities in opposite 
hemispheres, due to kinematical or geometrical effects, for example. R q is the hadronic 
branching fraction of the Z° to quarks q: 



Rq 



had 



They are taken to be the Standard Model values [21 



In ]T(J several so-called "hadronisation relations" based on approximate SU(2) 



symmetries between the different rj l q were proposed, such as rjtr ^77^. These extra 
constraints are necessary to solve the system of equations in the case of limited statistics, 
as at LEP. At a high luminosity e + e~ collider running at the Z° peak a solution can 
be found with fewer hadronisation relations, although w 77^ has to be kept. Note 
however that a systematic uncertainty of up to 2% should be assigned to this relation, 



following studies with QCD models frOfl . 



Although the relation i]^ ~ is fairly well motivated, it uncertainty would become 
a limiting factor in the determination of the branching ratios of the W boson and would 
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also make the measurement somewhat model dependent. We prefer to abandon the use 
of the relation by also taking into account the W determination. The W decays provide 
a separation of up and down quarks. At the Z° the main handle to separate and rfc 
is to tag an up or down event by a high x p pion in one jet and to find a charged kaon 
(indicating an up quark) or a neutral kaon (indicating a down quark) in the opposite jet. 
Given the dilution of the signal by decays of strange vector mesons and strange quark 
events, the discrimination power is only marginal. The resulting uncertainty in the if q 
limits the accuracy with which the CKM matrix elements can be determined. On the 
other hand, since the charged W boson can only decay into a restricted number of quark 
combinations, these decays discriminate between up and down quarks. For example, a 
charm quark can only be associated with a down quark, but not with an up quark. As a 
consequence, a combined fit of W and Z° decays for the r\ l and the W branching ratios 
allows one to find a solution without any assumptions about QCD symmetries. 



3.2 Experimental aspects of running at the Z° 

In addition to secondary vertex finding, the main selection requirements at the Z° are to 
identify hadrons. In order to have sufficient particle separation power with the dE/dx 
measurement and for the V° reconstruction of Kg and A identification, we require jets 
and tagging particles to be within the central part of the detector. Specifically we require 
that the polar angle of the thrust axis with respect to the beam direction, ^Thrust, satisfy 
I cos thrust I < 0.8 and that the tagging particles have a polar angle of the momentum, 9 P , 
within the range | cos# p | < 0.9. The detailed reconstruction criteria will depend on the 
performance of the tracking system. 

For the luminosities expected at TESLA at the Z°, we estimate from a JETSET 
simulation the number of single- and double-tagged events listed in Table ||. As discussed 
in |10| , |TT| and mentioned in the previous section, the tagging efficiencies r\ % q can be 
determined from these measurements with errors as given in Table |j. Also shown in 
Table f§ are the contributions from a ±2% systematic uncertainty in the hadronisation 
relation if A « rj^ . In omitting this relation but including data from W decays, the results 
are independent of any assumption about hadronisation and represent a considerable 
improvement over the precision reached at LEP. In general the correlations between the 
various elements are small. The most important correlations are between the rj l q of the 



same tagging particle type. For example, 77J and rf^ are almost fully anticorrelated. 



3.3 Systematic uncertainties 

In view of the unprecedented number of Z° decays available at GigaZ a detailed evaluation 
of the systematic uncertainties is unrealistic at this stage. We will just mention some 
potential distortions and suggest how to estimate them. From past experience we assume 
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particle 
type 


tagged 
hemispheres/1000 


7r± 


K± 


doub 
P(F) 


le-tagg 


ed even 
A(A) 


ts /1000 
charm tag 


bottom tag 




461395 


41890 


35442 


8667 


4442 


1962 


29183 


949 


K± 


222445 




10043 


4212 


2793 


1227 


16361 


546 


PIP) 


50187 






490 


542 


240 


3283 


110 


Kg 


29304 








204 


178 


2137 


69 


A(A) 


12747 










39 


831 


27 


charm tag 


587044 












124575 


51628 


bottom tag 


450683 














108898 



Table 3: Numbers of tagged event hemispheres and double-tagged events, scaled down by 
a factor of 10 -3 , for x p > 0.2 in an event sample of 2 x 10 9 Z° events from a JETSET 
simulation. 



that the huge amount of data will provide enough cross checks to keep all these sources 
of uncertainty under control. 

One crucial element of the analysis is charm and bottom quark tagging based on 
secondary vertices. As shown at LEP and SLD many contributions to uncertainties of 
their efficiencies and purities can be derived from data. 



As discussed in [[11] the major systematic uncertainties in the light quark sector in 
the LEP analysis are due to the efficiencies and purities of the hadron identification. At 
LEP these are estimated from relatively pure samples of particles from Kg — > 7r + 7r _ and 
D° — > K~7r + decays or photon conversions into an e + e~ pair, for example. Such cross 
checks can also be performed at GigaZ where one can expect the much higher statistics 
to lead to a sizeable improvement of the systematic uncertainty compared to LEP. 

Another major uncertainty in the LEP analysis comes from the hadronisation relations, 
briefly mentioned in Section \S.1[ These relations were needed to obtain a stable solution 
of the equation system. As discussed above, with the higher statistics at GigaZ and the 
use of W decays these hadronisation relations are no longer needed. 

Depending on the actual detector performance it may also be possible to use additional 
high x p particle types like 0(1020) mesons, which are very likely to originate from a strange 
quark, to further constrain the tagging efficiencies. 



4 Determination of the W branching ratios 



4.1 Formalism 

Neglecting experimental effects, the tagging efficiencies in the rest frame of the W and Z° 
bosons are almost identical, since the masses are so similar. Therefore, the rf determined 
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quantity 


value 


error 
fitting only Z° 

Vd ~ 


error 
fitting Z° and W 
no had. rel. 


vf 


0.209478 


0.000045 


0.002568 


0.000800 


vf 


0.208519 


0.000045 


0.002568 


0.001058 


vf 


0.129708 


0.000089 


0.000076 


0.000087 


vf 


0.028584 


0.000014 




0.000014 


vf 


0.000472 


0.000003 




0.000003 


'Id 


056490 


0.000276 


0.001558 


000475 


7u 


0.074605 


0.000437 


0.002761 


0.000819 


vr 


0.122413 


0.000089 


0.000594 


0.000175 


vr 


0.019890 


0.000011 




0.000011 


vr 


0.000247 


0.000002 




0.000002 


p 


0.014970 


0.000077 


0.000166 


0.000091 




0.024974 


0.000107 


0.000207 


0.000113 




0.019905 


0.000043 


0.000594 


0.000095 




0.003399 


0.000005 




0.000005 


p 

% 


0.000058 


0.000001 




0.000001 


K s 

Vd 


0.007233 


0.000089 


0.000445 


0.000146 


Kg 


0.005826 


0.000119 


0.000514 


0.000176 


Kg 

r)s 


0.019489 


0.000021 


0.000047 


0.000024 


K s 


0.002573 


0.000004 




0.000004 


Kg 


0.000029 


0.000001 




0.000001 




0.003048 


0.000052 


0.000183 


0.000069 




0.002893 


0.000072 


0.000228 


0.000091 




0.008498 


0.000015 


0.000007 


0.000014 




0.000858 


0.000002 




0.000002 




0.000014 


0.000001 




0.000001 



Table 4: Values of the and their expected precision for x p > 0.2. Column three gives 
the statistical uncertainties from a fit to 2 x 10 9 Z° decays assuming the hadronisation 
relation v(l i]^ . The fourth column shows the effect of a ±2% systematic uncertainty 
in the hadronisation relation. Column five gives the results of the combined fit to the Z° 
and W decays without invoking a hadronisation relation. 
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up 


charm 


down 


2375000 


126000 


strange 


124800 


2370000 


bottom 


200 


4000 



Table 5: Number of produced W boson decays into specific quark pairs for a total of five 
million accepted hadronically decaying W bosons. 

at the Z° peak can be used to measure flavour production in W decays in the W rest 
frame. To determine the branching ratios of the W bosons we use both single and double 
tags. A singly tagged W is where a tag of type k is found in just one of the jets. A doubly 
tagged W is a candidate where particle types i and j are tagged in each jet belonging to 
a W boson. This leads to the generic equations: 

N k = < ad -Ek fc ( 1 -E^) + ( 1 -E^]-^' ( 3 ) 

qq' I I 

Nq = iC-(l-0.5MS^ + #-fl M ' , (4) 

qq' 

where 5ij = 1 if i = j, and zero otherwise, which avoids double counting if identical 
tags are required, and where N^ d is the number of selected hadronically decaying W 
candidates. Here we neglected experimental effects, hemisphere correlations, background 
from non-W events, and assumed the right assignment of particles and jets to W bosons. 
In a real experiment none of the conditions holds exactly. The full equations which 
account for these complications are given in the Appendix. The number of 31 possible 
equations over-constrain the six unknown branching ratios, which can be obtained from 
a x 2 fit- 
As discussed in Section [3.1| , to minimise the uncertainties originating from the limited 
knowledge of the rj 1 , we perform a simultaneous fit for the r] l q and the W branching ratios 
using both Z° and W decays. 

4.2 Events with W Bosons at TESLA 

The high luminosity of TESLA leads to a substantial yield of W bosons. The main 
production processes are e + e~ — ► ez/ e W, with a cross section of some 5 pb at y/s = 
500 GeV, and W pair production of some 8 pb. In addition, W bosons are produced in 
top decays, but because of the lower cross section and some complications arising from the 
multi-jet environment we will not consider them. In total some seven million W bosons 
will be produced in a nominal TESLA year. 

We restrict this analysis to W bosons that are scattered into the central part of a future 
TESLA detector because of the experimental requirements of good hadron identification 
and efficient and pure heavy quark tagging. Within the polar angles 6*w of the W boson 
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\cos6w | < 0.8 some 700000 hadronically decaying W bosons can be retained for W 
pair production, while some 450000 single hadronically decaying W bosons are kept. We 
will base the following discussion on 5 million usable hadronic W decays, which could be 
collected in a few years of high energy data taking. Assuming the current knowledge of 
the CKM matrixfj, the expected yields of the different decay modes are listed in Table 

4.3 Experimental Procedure for W Bosons 

The basic strategy is to boost the decay products of W bosons into the W rest frame, 
search in each jet for a secondary vertex or for a light flavour tag. In the latter case we 
identify the type of the highest momentum particle in each hemisphere, and finally assign 
the probability rf obtained from Z° decays that such a tag stems from a certain primary 
quark flavour. To achieve this one has to first find the proper association of particles to 
W bosons, reconstruct the W boson energy and momentum and identify the tag. 

In case of single W production the association of particles to the W boson is 
unambiguous. This is also true for W pair production where one W decays hadronically 
and the other into a pair of leptons, denoted " semi-hadronic" . If both W bosons from a 
pair decay hadronically, the jet assignment is more difficult. After grouping the particles 
into four jets, three different jet pairings are possible. Our simulation studies using 
PYTHIA || show that at least for LEP energies this pairing is correct for about 85% of 
all W pairsQ. 

In W pair events the energies and momenta of the W bosons can be rather precisely 
reconstructed by fitting the observed momenta and energies of jets and leptons to obey 
energy and momentum conservation and to combine to have the W mass. This has been 
shown by the LEP experiments, which use total centre-of-mass energy as given by the 
accelerator and the fact that the W pairs are produced at rest at e + e~ colliders. Similarly 
in single W events the two jets can be constrained to have a mass identical to the mass 
of the W, M w . 

To identify the primary quark flavour, each jet is searched for a secondary vertex and 
classified either as bottom, charm or, if it does not have a secondary vertex, as a light 
flavour candidate. The efficiencies for the charm, bottom and light flavour tags are listed 
in Table |TJ. For the light flavour tags the tracks and clusters assigned to each W boson 
are then boosted along the reconstructed four momentum of the W candidate into its rest 
frame. The particles are grouped into two hemispheres with respect to the direction of 
the thrust axis and in each of these the particle with the highest momentum is identified. 

4 Reference 0| gives directly measured values and the ranges from a unitarity constraint CKM matrix. 
For our purposes we adjusted these values such that the individual decay modes add up to exactly five 
million W bosons. We also assumed the decay into up quarks to be exactly as frequent as the one to 
charm quarks. For the purpose of this analysis these adjustments are unimportant. 

5 No detailed study exists yet for energies of 500 GeV. We therefore assume for the following evaluation 
the purity observed at LEP. 
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particle 
type 


tagged 
hemispheres 


_|_ 


t rA- 


doi 
P(P) 


rble-taj 
K s 


$ged ev( 
A(A) 


arts 
charm tag 


bottom tag 




829270 


116657 


84621 


23297 


9287 


4089 


240508 


7639 


K± 


349183 




17061 


8556 


4035 


1714 


190195 


6300 


PIP) 


87307 






1116 


997 


429 


33204 


1077 


K° s 


41863 








235 


198 


29319 


980 


A(A) 


18241 










41 


12775 


426 


charm tag 


1124876 












62888 


4517 


bottom tag 


36701 














81 



Table 6: Numbers of single and double tagged candidates obtained from a sample of five 
million W decays from a PYTHIA simulation. 



These particles are retained if they are in the geometrical acceptance range, fulfil the 
criteria to identify their species and have a scaled momentum in the W boson rest frame 
x p = 2p/Mw > 0.2. The expected numbers of single and double tags in five million W 
decays are listed in Table ||. 

To apply the Z° tagging efficiencies to the W decays, several corrections have to be 
applied, even though the selection and reconstruction resembles closely the one used to 
determine the rf at the Z°. 

• The mass of the W is some 10 GeV below the Z° mass. For particle tags this 
leads to a small difference of the energy spectra due to QCD scaling violations. 
The resulting differences in the tagging efficiencies at these two energy scales were 
estimated with e + e~ continuum events generated with JETSET at centre-of-mass 
energies equivalent to the Z° and the W masses. The fractional changes of the rf q 
relative to those at the Z° are shown in Figure [I]. For x p > 0.2 the r] q were found 
to be 1.7% higher at the W mass. These scaling violations are rather independent 
of the primary quark flavour. Some dependence on the tagging particle is observed. 
Note that the latter can be rather precisely determined from data by summing over 
all flavours. 

Furthermore, differences in the detection efficiencies have to be accounted for: 

• The high x p particles and those from the secondary vertices from W and Z° decays 
have different momentum and polar angle distributions in the lab system. This can 
be inferred from Figures ^| and |3| where the jet energies in single W and W pair 
events are shown. There is a broad spectrum of jet energies, mostly larger than the 
jet energy of ~ 45.5 GeV in Z° decays. This potentially leads to different efficiencies 
and purities of secondary vertex finding and, for the highest x p particle, momentum 
resolution and dE/dx efficiencies. 



11 





Hadronic branching ratios R qq > 


CKM matrix elements \V qq >\ 




up 


charm 


up 


charm 


down 


0.4750±0.0027 


0.0252 ±0.0016 


0.9747±0.0028 


0.2245 ±0.0072 


strange 


0.0250±0.0027 


0.4740 ±0.0017 


0.2234±0.0124 


0.9737 ±0.0017 


bottom 


0.00004±0. 00010 


0.00080±0. 00005 


0.0089±0.0114 


0.0400 ±0.0011 



Table 7: Expected precision on the hadronic branching ratios of W bosons and CKM 
matrix elements. Note that R u ^ is constrained by R u ^ = 1 — J2 qq > Rqq'i where the sum 
ranges over all other hadronic branching ratios. 



• In addition possible misassignments of particles to W bosons and distortions in the 
reconstruction of the momentum of the W bosons have to be taken into account. 



The necessary corrections can be estimated using detector simulations. In many cases 
significant cross checks with data can be performed. In applying the r] q to fully 
hadronically decaying W pairs, there could emerge in principle a further distortion due 
to Bose-Einstein or colour reconnection effects. However, as already known from the 
discussions on the colour reconnection for W pairs at LEP, this affects in particular low 
momentum particles, whereas the leading particles are rather undistorted 

In addition background from non-W production processes have to be considered. The 
major backgrounds in the W pair sample with fully hadronic decays are due to e + e~ 
continuum quark production with two hard gluons, Z° pair production and top pair 
production. Early studies at 500 GeV indicate that these backgrounds can be kept 



below the 5% level. Background is significantly smaller for semi-hadronic events. 



5 W Branching Ratios and CKM Matrix Elements 
5.1 Statistical precision 

Following the procedure outlined above we estimate the precisions for the hadronic 
branching ratios, R qq >, using a \ 2 fit to the observed tagged yields in W and Z° decays, 
in W and Z° decays. In the following we assume that the W decays only into the known 
leptons and quarks, except the top quark which is inaccessible because of its mass. If 
additional decays contribute they probably will either modify the true R qq > or even imply 
that Y2-R-qq' ^ !• As default we assume the Standard Model case in this analysis and 
constrain i? u b = 1 — J2 qq > Rqq'i where the sum ranges over all other hadronic branching 
ratios. The result of this fit is listed in Table |7j. We just note that by omitting such a 
constraint we can determine Y2 qq ' Rqq' with a statistical precision of ~ 5 x 10~ 4 . 

Without assuming unitarity, the CKM matrix elements can be obtained from the 
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partial widths of the W bosons, in a way similar to what has been used in ||: 

r(w^ m ') = c ^^ \v„,\\ (5) 

where Gf is the Fermi constant and 



Cqod = 3 {I + <M^1 + i.409 (Z&ST) 2 - 12.77 (^>) ' > (6) 

expresses the QCD radiative corrections. The partial width is related to the hadronic 
branching fractions by: 

T(W - qq') = T w B h R qq , (7) 

with T w the total W boson width and B h the inclusive branching ration of the W into 
hadrons. 

The CKM matrix elements can thereby be determined using the measured values of the 
basic quantities. However, such a measurement may finally be limited by the knowledge 
of Tw- Given by the very precise measurement of Gf in \x decays and lepton universality 
measurements in r decays ||25|| , the latter can be substituted by the measurement of the 
electronic or muonic partial width of the W leading to 

F(W -> qq') _ T w B h R qq , _ 2 

F?W 7 _ ~P D — ^QCD ■ \Vqq'\ , \°) 

where Bu tfl \ is the electronic or muonic decay fraction of the W. Thus, the CKM matrix 
elements can be related to the measured fraction of W — > qq' decays using the QCD 
correction factor and the well-measurable ratio B^/B^^ 



\ V qq'\ = 1T^- R qq' ■ ( 9 ) 

<^QCD -D(e,At) 



Currently the ratio B^/Bu^) is known experimentally from LEP data to 1%, but the 
error should significantly decrease to (9(0.05%) due to the higher yield of W bosons at 
TESLA. The QCD radiative correction is calculated to third order in as', however, the 
strong coupling is currently only known to a few percent. Given two billion Z° events the 
uncertainty could be rather significantly reduced. In fact Cqcd should be almost identical 
in Z° and W decays. Both these error contributions are smaller than the uncertainties in 
Rqq 1 - We obtain the matrix elements given in Table and their correlations are listed in 
Table |^. We do not show the negligible correlations between the \V qq i\ and the rf resulting 
from the combined fit. The results are consistent with the input values. 

At this stage we have completely neglected any possible theoretical corrections and 
uncertainties due to mass effects, electroweak contributions etc. 
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1 t r 1 

1 Kid 1 


I t r l 

IKs| 


1 t r 1 

IKbl 


It/ 1 

IKdl 


It/" 1 

IKsl 


It/" 1 

IKbl 








1.000 


-0.983 


0.087 


-0.590 


0.566 


-0.014 




Ks 






1 C\C\C\ 

1.000 


-0.174 


0.581 


-o.60o 


0.006 




Kb 








1.000 


-0.117 


0.139 


-0.107 




Kd 










1.000 


-0.964 


0.014 




Ks 












1.000 


-0.014 




Kb 














1.000 



Table 8: Correlations between CKM matrix elements. 



5.2 Systematic Uncertainties 



For the envisaged precision and the huge amount of data, a detailed analysis of the 
potential systematic uncertainties is impossible at this stage. Here we discuss some major 
uncertainties and indicate ways to estimate their importance. In all cases we believe that 
they can be kept under control. 



Errors and correlations for 77* from Z°. The uncertainties of the efficiencies 
r] l q from the billion Z° bosons produced at a high luminosity GigaZ propagate into 
an uncertainty in the branching ratios of the W. This is taken into account by the 
result of the combined fit. 

The principal problem will be to keep the systematic uncertainties of the rf q under 
control to this high level of precision. Some potential uncertainties were mentioned 
in Section |3~3"|. Ideally the data at GigaZ and at high energies would be collected 



with the same detector. In this case many of the potential uncertainties are the 
same for the Z° and the W boson analyses and can be neglected. Otherwise the 
relative uncertainties between two different detectors in flavour tagging have to be 
taken into account. Given the large statistics of data to cross check the detector 
performance, it appears possible to maintain a sufficiently high precision. 

W boson reconstruction. The accuracy of the measurements of the branching 
ratios, R gq >, depends on how reliably particles and jets can be assigned to W 
bosons and how well this is understood. Related is the correctness of the W boson 
reconstruction from a kinematic fit. These affect the determination of the x p values 
of the tagging particle in the W rest frame. 

In part potential misassignments can be estimated from other processes at high 
energies, like e + e~ — > Z°Z°. Our simulation studies suggest that the corresponding 
error will be unimportant. 



Differences in tagging efficiencies. As discussed in Section |4.3| , several subtleties 
have to be taken into account in translating the tagging efficiencies at the Z° to those 
for the W bosons. The differences can be minimised by applying rather similar 
selection requirements in the two cases. 

Here we list potential corrections: 
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— relative capability for charm, bottom and light flavour tagging (including 
particle identification capability) in W pair production at high energies 
compared to the Z° data, 

— QCD scaling violations which slightly change the rj l , 

— effects of the different angular acceptance and kinematical and geometrical 
correlations. 

None of these uncertainties appears to impose a substantial problem. Systematic 
uncertainties can be derived from data. 

• Backgrounds. Most of the overall yield from potential background processes will 
be measured to very high precision and their contribution to the sample of W bosons 
can be well determined. Somewhat more uncertain is their contribution to the light 
flavour tags. This requires a knowledge of the particle content. For example, the 
number of e + e~ — ► qq continuum events with two hard gluons will probably be 
quite well known. However, to account for their contribution to the W sample, 
the leading particle distributions in gluon jets also needs to be understood. A 
possible way to determine their impact is to study the particle composition in these 
processes for kinematic properties that cannot be confused with W pair production 
and extrapolate it to the relevant kinematic configuration. 



5.3 Some possible improvements 

• Particle separation. A better detector performance to separate the hadron species 
will lead to a better purity of the individual light flavours. 



Tags including charge. As already discussed in |10[, including the charge sign of 



the tags allows one to improve the flavour purity of the 77* from Z decays. Charge 
dependent rf q are even more interesting for both W production processes because of 
the significant charge asymmetry in cos#. 

Polarisation. W pair production at TESLA will lead to polarised W bosons, the 
degree of polarisation depending on the scattering angle, leading down-type jets 
to be in general more energetic than up-type jets. This effect can be seen in the 
accumulation of jets with high and low energies in Figure |2|, for example. One 
may use this property to statistically separate between up- and down-type quarks. 
Ambiguities from the high x p tagging between charged kaons in up and strange 
quark events, for example, can thereby be reduced. Such separation could become 
even more powerful if in addition the charge of the tagged particle is used. 
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6 Top quark decays 



As yet the CKM matrix elements have only been determined indirectly from processes 
such as b — > S7 and B°B° mixing. In both cases the transition involving top quarks occur 
at the one loop level. A first direct measurement of the fraction of top decays into bottom 
quarks can be found in |2fJ. The high statistics and clean environment at a high energy 
e + e~ collider allow more comprehensive and precise direct measurements. 



6.1 Formalism 



To extract the relevant branching ratios of top quarks into the various species of down - 
type quarks one has to tag just the jet that is produced together with the W. An equation 
system can be constructed relating the number of events, A^, identified with a tag k, out 
of a number of top decays, iVtop, 

N k = N top V k q B(t^qW) (10) 

g=d,s,b 

which can be easily solved for the branching ratios B(t — > qW). In the following we 
assume the absence of any non-standard decay. As was discussed for W decays, such an 
exotic contribution may become apparent by measuring the branching ratios with high 
precision. 



6.2 Top Yields at TESLA 



With the luminosity given in Section |^, some 400,000 top quarks will be produced each 
year. To determine the transition t — > qW it is least ambiguous to demand leptonic W 
decays, thus avoiding any confusion in jet assignments. With such a requirement, some 
120000 top quarks per year will be kept. However, it may also be possible to use at least a 
part of the fully hadronic top decays, providing a significantly larger sample. We will base 
our discussion on one million well-tagged and reconstructed top quarks in the W — > lu 
mode. The expected yields of the various decays are listed in Table || 



6.3 CKM Matrix Elements from Top Quark Decays 



The branching ratios of the top quark are related to the CKM matrix elements by including 
higher order corrections summarised in P7[ : 

r(t -> qW) ~ \V tq \ 2 x 1.42 GeV . (11) 
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|V*| 


number of top decays 


down 


0.006 


40 


strange 


0.04 


1600 


bottom 


0.999 


998360 



Table 9: Number of expected top quark decays into down type quarks. Here we assume 
one million top decays. Also shown is the CKM matrix element assumed. As for the 
estimate on W decays we fixed the CKM matrix elements to be consistent with unitarity. 



It should be noted that, in contrast to all other known quarks, the top quark is expected 
to decay before it hadronises This renders the determination of the CKM matrix 

elements independent of uncertainties inherent in hadron physics. 

The partial width r(t — > qW) can be related to the branching ratios B(t — > qW) by 

r(t->gW) = T top B(t^qW) (12) 



and thus 



'^l 2 = B < 4 - TOT ■ < 13 > 



Whereas the determination of the branching ratios is straight-forward and will be 
detailed in the next paragraphs, the measurement of the decay width of the top quark to 
high precision is less obvious. This is not possible via secondary vertex tagging because 
of the expected lifetime of O(10~ 25 sec). Instead other means have been suggested |29|| , 



the most promising one seeming to be single top production e + e~ — > ez/ e tb P0| . For 
such a measurement the background from top pair production has to be rejected. Scaling 
the precision estimated in |3U] to the luminosities foreseen at TESLA, r top should be 
measurable to ~ 1%. 

Assuming unitarity and that the top quark couples only to down, strange and bottom 
quarks, the r top should agree with the theoretical expectation and thus 

\V tq \ 2 = Bit -+ qW) (14) 



In general the determination of the branching ratios follows the procedure for the W 
bosons. 



• Identification of the W boson decay products and consideration of the remaining jet 
in a top decay. This is rather straight-forward for leptonic W decays. 

• Application of the tagging algorithm to this jet. The light flavour tag should have a 
very high purity even at the cost of a lower efficiency because of the overwhelming 
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7T* 


K ± 


P(p) 


Kg 


A(A) 


bottom tag 


Number of tags 


1000 


620 


125 


90 


50 


499180 



Table 10: Numbers of expected top quark decays tagged by a light quark tag combined 
with a high x p particle or a heavy quark tag. 





Branching ratios 


CKM matrix elements 


t -> d 


(8±52) x 10~ 5 


0.0060 ±0.026±0.00003 


t -> s 


0.0015± 0.0005 


0.0400 ±0.006±0.0002 


t->b 


0.99840±0. 00028 


0.999200±0.000008±0.005 



Table 11: Expected precision on the CKM matrix elements from top decays assuming 
top only decays into d, s, and b quarks. The second error of the CKM matrix elements 
assumes r top to be known to 1%. 



fraction of top decays into bottom quarks. In the following we assume e u d s = 0.5 and 
eb = 0.01. In the absence of FCNC decays, no charm quarks are directly produced. 
This allows one to reject bottom quarks by a very tight selection against secondary 
vertices, helping to achieve the preferred high purity of light flavours. 

• Different from the situation with W decays, the mass of the colour neutral system 
in which hadronisation takes place depends on the details of the decays of both 
the top and anti-top quarks in the same event. This makes the assignment of the 
if q somewhat more complicated. Whereas the Z° and the W bosons are colourless 
objects and therefore hadronisation proceeds in their respective rest frames, the top 
quark itself is coloured. Its colour is neutralised by the recoiling anti-top quark. 
In general the hadronisation should evolve in the qq' system given by the decays 
(t — > gW + ,t — > q'W~). This requires that the mass M qq i of the qq' system be 
reconstructed. QCD scaling violations between the rf at the Z° and those at M qq i 
have to be taken into account. A JETSET simulation for M qq > is shown in Figure |] 
and the sizes of the necessary QCD corrections are given in Figure [l| 

If at least one of the top quarks decays fully hadronically, distortions may arise from 
colour reconnection between quarks from the W decay and those directly associated 
to the top. As mentioned previously, they are expected to be of no importance for 
our method. The experimental task, however, is complicated for hadronic W decays 
by the need to resolve the two jets that directly couple to the top quarks. This may 
require further selections. 

• Solution of the equation system (Eq. |i~0|) . 



The numbers of expected tags are listed in Table 10 



The expected precision of the 
CKM matrix elements involving top quarks are listed in Table [IT]. Their correlation 
matrix is given in Table [T2 . 
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\Vt d \ 


\Vt B \ 


|Kb| 


IKdl 


1. 


-0.832 


0.728 


\Vts\ 




1. 


-0.986 


\Vtb\ 






1. 



Table 12: Correlation matrix between the top CKM matrix elements. 



These results imply that one million top decays allow one to determine the branching 
ratio B{\,—* bW) to be different from unity with high significance. This implies that 
the sum of Bit — > sW) and B(t — > dW) is larger than zero. The light flavour tagging 
efficiencies even allow one even to distinguish down and strange quark contributions. As 
a result the CKM matrix element |Vt s | can be determined directly and with reasonable 
precision. Despite low statistics and still formidable bottom background in the light 
quark samples, a first upper limit from direct measurements on the t — > d transition can 
be derived. At 95% confidence, 

B(t -> dW) < 8 x 10" 4 . 

To translate these into reasonably precise |Vt g | elements requires that r top be measured 
with higher precision than what has been suggested so far. If this can be achieved one 
can set a limit of 

\V td \ < 0.029 . 



7 Conclusions 



The large samples of W bosons and top quarks at a future e + e~ collider of high energy 
and luminosity provides the means to directly determine the fundamental CKM matrix 
elements. The current values are based on hadron decays and deep inelastic scattering 
results and invoke QCD symmetries. We propose a complementary method which is free 
of any assumptions on QCD modelling at a low mass scale. In addition the resulting 
precision is at least competitive to the ones from hadronic decays. In this paper we have 
only discussed the experimental feasibility. The potential experimental precision may 
have to be complemented by further theoretical scrutiny of higher order effects. 



Compared to the current results for individual measurements as summarised in |22| 



there are improvements on all elements in the charm sector. The precision is better than 
what is anticipated from other measurements in the future. Even measurements of |V^b| 
at B factories hardly reach the potential TESLA precision because of the anticipated 
theoretical uncertainty. The results for transitions involving the up quark are a factor 
three to five worse than those from hadron decays. However, even for these elements the 
measurement in W boson decays may become an interesting complement. 

The direct determination of CKM matrix elements involving top quarks may be 
pioneered at TESLA. So far only |Vtb| has been determined directly using top decays 
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nil T*T/3'n"l" n n r^t^y^ - o l Y\~t"\r 
LyLlllclllj LlllL-Cl LcLlllLV 












± 0.0008 




±0.0028 




v us 




i o nnoQ 
ztU.UUzo 




±U.Ulz4 




Kb 




~ ±0.008 


±0.0004 Q 


±0.011 




v cd 




±U.UlD 




±U.UU ( z 




v cs 




±0.16 




±0.0017 




Kb 




±0.0019 


±0.0012 [0] 


±0.0011 




Kd 




|Kd|/|K| <0.24 


<0.016 


±0.026 ±0.00003 




K 




±0.006 ±0.0002 




Kb 




(±0.29,-0.12) |2C| 


~ ±0.05 |27| 


±0.000008 ±0.005 



Table 13: Current (mostly from |[22|| ) and expected precision of measurements of the CKM 
matrix elements. Only direct measurements and those that do not rely on the knowledge 
of other CKM matrix elements are used. Prospective measurements are listed for TESLA, 
B factories [31j and the LHC [27|. For the TESLA limits on |Kj| the second error 



e 1 e 



is due to an uncertainty of r top of 1%. 



at the Tevatron [26]. The measurement of the branching ratio B(t — > bW) may be 



improved at the Tevatron Run 2 and particularly at the LHC. For example, [£7[ assumes 
a 0.2% statistical error for one year of low luminosity running at LHC (10 fb _1 ). The 
systematic error is not yet evaluated. Translating this into |Kb| is more complicated 
and a direct measurement via single top quark production is deemed "challenging". No 
direct measurements of |Kd| and |K| are possible at the LHC. At TESLA not only the 
dominant decay branching ratio B(t — > bW) can be measured with similar or better 
precision compared to the LHC, but also the top transition to light quarks can be probed. 
The element |K| can be determined with a significance of about six standard deviations 
and, if r top is known to sufficient precision, a significant upper limit on |Kd| can be set. 
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A Formulae for Single and Double Tags in W Decays 



Equations ^] and § show only the general structure of the equation system to solve for the 
branching ratios of the W bosons. As pointed out, these equations become more involved 
if experimental distortions are taken into account, such as kinematical and geometrical 
correlations pw between the two jets, and the probability n w ^j et of correctly assigning a 
jet to a W boson. 
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For A^ d denoting the number of accepted hadronically decaying W bosons, the 
number of double tags of types i and j with a correct assignment of a jet to a W boson 
is given by 

^correct = ^had . . (1 _ . 5 ^)PW J>X' + rf^,) R qq , (15) 

qq' 

where the sum ranges over all six quark pairs possibly produced in W decays: qq' = (ud, 
us, ub, cd, cs, cb). In W pair events where the two W bosons decay like Wi — > q± q~\ and 
W 2 — > g 2 g 2 , those with misassigned jets can be approximated by 



^assigned = < ad (l - n w ^ jet )(0.5 - 0.25*y) 



9191 

.9191 



^^92+^92)^9292 
.9292 

(16) 



In this case an incoherent sum over the branching ratios is assumed. Correlations 
between wrong combinations slightly modify the equation. As an example assume a 
true W + W~ — > (cs)(ud) event. If, for example, the c quark is combined with the d quark, 
the other combination is fixed as being su. Such a correlation is not accounted for in the 
above equation. However, it is rather straight-forward to include it. Such a modification 
is important only if all four quarks are tagged, which is rather unlikely. 

In addition one has to estimate the double tags from background, mainly from 
continuum QCD events with two hard gluons and from Z° pair events. The amount 
and particle content of these events must be estimated from simulation. 

The single tag rate is given by: 

N t = N w { 

rii [IC d -B{W -> ud)+/C s -5(W -> us)+/C b - J B(W -> ub)] 
+ /C u [n d -B(W -> ud)+ril -B(W -> us)+r£ -B{W -> ub)] (17) 
+ 4 [lC d -B(W -> cd) +/C S -5(W -> cs) +/C b -5(W -> cb) ] 
+ /C c [r] d -S(W -> cd) +^ -S(W -> cs) +r] l h -B(W -> cb) ] } 

where /C 9 is the probability that jet from quark g does not lead to any tag: 

JC q = [1 - p w + (1 + Pw)^<]. (18) 

9 
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Figure 1: Fractional difference (rf q (E cm ) — rf q (M z o)) / 'r] l q (M z o) for various if q as a function 
of the centre-of-mass energy, E cm , from the JETSET Monte Carlo. 
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Figure 2: Energies of jets in W pair events generated with the PYTHIA Monte Carlo at 
a centre-of-mass energy of 500 GeV. 
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Jet Energy in Single W Events (GeV) 



Figure 3: Energies of jets in single W events generated with the PYTHIA Monte Carlo 
at a centre-of-mass energy of 500 GeV. 
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qq' Invariant Mass (GeV) 



Figure 4: Invariant mass of the qq' system from the decays t — > gW + and t — > q'W 
generated with the JETSET Monte Carlo at a centre-of-mass energy of 500 GeV. 
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